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Abstract. 

QQ , We study black hole solutions in general relativity coupled to a unit timelike vector 

^S| ' field dubbed the "aether". To be causally isolated a black hole interior must trap 

matter fields as well as all aether and metric modes. The theory possesses spin-0, spin- 
1, and spin-2 modes whose speeds depend on four coupling coefiicients. We find that 
QQ ' the full three-parameter family of local spherically symmetric static solutions is always 

OO . regular at a metric horizon, but only a two-parameter subset is regular at a spin- 

^-p I horizon. Asymptotic flatness imposes another condition, leaving a one-parameter 

f^ ■ family of regular black holes. These solutions are compared to the Schwarzschild 

\^ , solution using numerical integration for a special class of coupling coefiicients. They 

^^ ' are very close to Schwarzschild outside the horizon for a wide range of couplings, and 

O . have a spacelike singularity inside, but diflFer inside quantitatively. Some quantities 

constructed from the metric and aether oscillate in the interior as the singularity is 
approached. The aether is at rest at spatial infinity and flows into the black hole, but 
differs significantly from the the 4-velocity of freely-falling geodesies. 
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1. Introduction 

Einstein- Aether theory, nicknamed "ae-theory", consists of general relativity (GR) 
coupled to a dynamical, unit timelike vector field u"" called the "aether". The aether 
is somewhat like a nonlinear sigma model field with a hyperbolic target space. But, 
unlike for a sigma model, the hyperboloid is dynamically determined at each point by 
the spacetime metric and the aether couples to the metric via covariant derivatives. 
The aether dynamics thus modifies the metric dynamics, and the coupled system differs 
significantly from the familiar examples of GR coupled to matter fields. 

Properties of ae-theory have been extensively studied for the past few years, and 
many observational signatures have been worked out (see [1] for a recent review). The 
primary motivation for studying ae-theory is to explore the gravitational consequences 
of local Lorentz symmetry violation. The aether itself breaks local boost invariance at 
each point, while preserving rotational symmetry in a preferred frame. There are four 
free coupling coefficients in the general ae-theory Lagrangian (besides the gravitational 
constant), thus providing a somewhat general class of boost violating model^il. There 
is a common range of the four coupling coefficients for which simultaneously all post- 
Newtonian parameters agree with those of GR [H El E], where all perturbations are 
stable |71 [8] and have positive energy [H [9l [10], gravitational and aether vacuum 
Cerenkov radiation [[llj is absent [6j , the gravitational constant in the Friedman equation 
agrees with Newton's constant [12], and gravitational radiation damping of binary pulsar 
orbits agrees with the rate in GR for weak fields at lowest post-Newtonian order ^U\ . 

In spherical symmetry the aether has one degree of freedom, a radial tilting. 
Restricting to the time-independent case, we found in [13] that there is a three-parameter 
family of local solutions, of which a two-parameter sub-family is asymptotically fiat. In 
[T3] we focused on the one-parameter sub-family of these in which the aether is aligned 
with the static Killing vector. These "static aether" solutions were found analytically, 
up to the inversion of a transcendental function. The solutions are characterized by 
their total mass, and it was shown in [13] that the solution outside a spherical star is in 
this family. 

In this paper we focus on the time-independent, spherically symmetric black hole 
solutions of ae-theory. This is the first step in comparing ae-theory with astrophysical 
observations of black holes. It is also of mathematical interest as an example of unusual 
black hole behavior with non-linear self-gravitating fields. 

The Lagrangian for ae-theory depends on four dimensionless coupling coefficients 
Ci,2,3,4- The static aether solutions referred to above depend only on the single 
combination ci + C4, but this is not the case for the black hole solutions. Moreover, 
unlike those solutions, the black hole solutions cannot (as far as we know) be obtained 
analytically, so all of the results in this paper are numerically obtained. We do not make 
an exhaustive study here for all values of the Cj, but rather just attempt to determine 
the generic behavior of the black hole solutions. 

I For a much broader class of models with Lorentz violating gravity see [21 [3] . 
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We begin this paper in section H] with a discussion of the definition of a black 
hole in ae-theory, and the conditions for the existence of regular black hole solutions. 
The qualitative reasons for the existence of a one parameter family of such solutions 
are explained. In preparation for the subsequent detailed analysis the action, field 
equations, and field redefinition properties of the theory are reviewed in section [3l This 
is followed in section H] by a demonstration using the power series solution of the field 
equations about a metric horizon that such horizons are generically regular, i.e. there 
is a three parameter family of solutions in the neighborhood of such a horizon, just 
as in the neighborhood of any generic point. We then show that spin-0 horizons are 
generically singular, but a two-parameter family of solutions is regular and the additional 
condition of asymptotic fiatness further reduces this to a one-parameter family. In 
section |5] we study the properties of typical examples of this one parameter family of 
black holes, imposing regularity by a power series expansion about the spin-0 horizon 
and numerically integrating out to infinity, determining the asymptotically fiat solutions 
by tuning the data at the horizon. These black holes are rather similar to Schwarzschild 
outside the horizon, and like Schwarzschild they have a spacelike curvature singularity 
inside at (or very near) zero radius. Unlike the static aether solutions of [13], m° is not 
aligned with the static Killing field in these black hole solutions. The aether fiows into 
the black hole, but differs significantly from the 4-velocity of freely-falling geodesies at 
rest at infinity. We also note that some functions constructed from the metric and aether 
exhibit oscillatory behavior in the interior as they approach the singularity. Section [6] 
concludes the paper with a brief discussion of questions for further work. 

2. General properties of black holes in ae-theory 

The relevant notion of a black hole in ae-theory is not immediately clear. To trap 
matter influences a black hole must have a horizon with respect to the causal structure 
of gab, the metric to which matter couples universally (or almost universally, according 
to observations). We call this a "metric horizon". This is not the only relevant notion 
of causality however. For general values of the coupling coefficients q, ae-theory has 
multiple characteristic hypersurfaces. In particular, perturbing around flat spacetime 
it was found in [7] that there are spin-2, spin-1, and spin-0 wave modes, with squared 
speeds relative to the aether given by 

spin-2 1/(1 -C13) 

spin-1 (ci - |ci + ^cl)/cu{l - C13) (1) 

spin-0 Ci23(2 - cu)/cu{l - Cis){2 + C13 + 3c2) 

where C13 = ci + C3, etc. The speeds Si are generally different from each other and 
from the metric speed of light 1. Only in the special case where C4 = 0, C3 = — Ci, 
and C2 = ci/(l — 2ci) do all the modes propagate at the same speed. In general, the 
characteristic surfaces for a mode of speed s, are null with respect to the effective metric 
Vab + (s^ — l)MaM6; whcrc r]ab is the flat metric and m^ the constant background aether. 
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In the nonlinear case characteristics can be defined as hypersurfaces across which the 
field equations admit a discontinuity in first derivatives [H] . For this paper we presume 
that these characteristics define the relevant notion of causal domain of dependence for 
the ae-theory field equations. This seems quite plausible, although no rigorous study 
has been attempted. The characteristic hypersurfaces are determined by the highest 
derivative terms in the field equations, so can also be identified by examination of high 
frequency solutions to the linearized equations about a given background [H]. For 
such solutions the gradients in the background are irrelevant, so we can infer that the 
nonlinear characteristics are null surfaces of the effective metrics, 

9ab = 9ab + (S- - l)UaUb. (2) 

We refer to the horizons associated with these metrics as the spin-0, spin-1, and spin-2 
horizons. If a black hole is to be a region that traps all possible causal infiuences, it must 
be bounded by a horizon corresponding to the fastest speed. The coupling coefficients 
Cj determine which speed is the fastest. 

A horizon is potentially a location where a solution to the field equation can become 
singular. This is because at a characteristic surface the coefficient of a second derivative 
term usually present in the equation vanishes. As such a surface is approached, the 
smallness of that coefficient may generically produce a solution in which some second 
derivative grows without bound, leading to singular behavior. This does not occur at 
spin-1 and spin-2 horizons in spherically symmetric solutions to ae-theory, presumably 
since there are no spherically symmetric spin-1 or spin-2 modes. However there is a 
spherical spin-0 mode, and we find that spin-0 horizons are generically singular. 

The requirement that the spin-0 horizon be regular reduces the three parameter 
family of local static, spherically symmetric solutions to a two parameter family, which 
reduces to a one parameter family when asymptotic fiatness is imposed. Hence there is 
just a one parameter family of regular static, spherically symmetric black hole solutions 
in ae-theory, just as in GR. Unlike outside a star, the aether in these solutions is not 
aligned with the Killing vector but rather fiows into the black hole. We conjecture that 
the collapse of a spherical star to form a black hole is nonsingular, and therefore must 
be accompanied by a burst of spin-0 aether radiation as the aether (and metric) exterior 
adjusts. 

The fact that the aether is not aligned with the Killing vector in a static black 
hole solution is no accident. It cannot be so aligned, since it is everywhere a timelike 
unit vector and the Killing vector is null on the horizon. Thus at a regular horizon the 
aether must be "infalling" , although it is nevertheless invariant under the Killing fiow. 
The static aether solution found in |13] has aether aligned with the Killing vector, and 
can be thought of as an extremal black hole with a singular horizon on which the aether 
becomes infinitely stretched. 

While the aether can be regular at a generic point on the horizon, it cannot smoothly 
extend to the bifurcation sphere B, i.e the fixed point set of the Killing fiow at the 
intersection of the past and future horizons [1]. The Killing fiow acts as a Lorentz 



Black Holes in Einstein- Aether Theory 5 

boost in the tangent space of any point on B, so it is impossible for the aether to be 
invariant under the flow there. This implies that the aether must blow up, becoming an 
infinite null vector as B is approached. This in turn raises the concern that there may 
be no regular metric horizon, since regularity on a future horizon is typically linked to 
regularity at B. Indeed Racz and Wald |15j have estabhshed, independent of any field 
equations, conditions under which a stationary spacetime with regular Killing horizon 
can be extended to a spacetime with a regular bifurcation surface, and conditions under 
which matter fields invariant under the Killing symmetry can also be extended. In 
spherical symmetry these conditions are satisfied for the metric, but the aether vector 
field breaks the required time reflection symmetry so it need not be regular at the 
bifurcation surface (although all scalar invariants must be, as must the aether stress 
tensor if the field equations hold) . 

Another potential obstruction to the existence of regular ae-theory black hole 
horizons arises from the form of the aether stress tensor. At a regular stationary 
metric horizon the Raychaudhuri equation for the horizon congruence with with null 
generator k°- implies that Rahk'^k^ must vanish, hence the Einstein equation implies that 
the matter stress tensor component Tabk°'k^ must also vanish. With common matter 
fields, e.g., scalar fields. Maxwell or Yang-Mills fields, and nonlinear sigma model fields, 
it is easy to show from examination of the form of the stress tensors that this condition is 
automatically satisfied locally for any field invariant under the Killing flow, independent 
of field equations. This property does not seem to hold kinematically for the aether stress 
tensor, but since we find a full three-parameter family of regular metric horizons, it is 
evidently imposed by the field equations. 

The fact that Tabk°'k^ does not vanish kinematically might appear to contradict 
the following general argument. For a Killing horizon with non-zero surface gravity 
it is not necessary to examine the form of the stress tensor to arrive at the inference 
that the horizon component of a matter stress tensor vanishes. If x° is the horizon- 
generating Killing vector, then the vanishing of the scalar TahX"")^ o^i fhe horizon is 
guaranteed by the facts that (i) it is invariant along the flow, (ii) x" vanishes at the 
bifurcation surface, and (iii) Tab is regular at the bifurcation surface (as guaranteed by 
the Racz- Wald extension theorem). But this argument too seems to fail for the aether 
stress tensor, because (as noted above) there is no purely kinematic way to argue that 
it (and therefore its stress-tensor) is regular at the bifurcation surface. So, again, the 
field equations seem to play an essential role in ensuring the existence of regular metric 
horizons. 

3. Einstein- Aether theory 

The action for Einstein-Aether theory is the most general diffeomorphism invariant 
functional of the spacetime metric Qab and aether field u"" involving no more than two 
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derivatives, 

where 

L = -R- K'^'mn'^aU'^VbU'' - \{gabU''u'> - 1). (4) 

Here R is the Ricci scalar, K°'^rnn is defined as 

where the q are dimensionless constants, and A is a Lagrange muhipher enforcing the 
unit timehke constraint. This constraint restricts variations of the aether to be spacehke, 
hence ghosts need not arise. A term of the form RabU°'u^ is not exphcitly included as it 
is proportional to the difference of the C2 and C3 terms in ([3]) via integration by parts. 

The metric signature is (H ) and the units are chosen so that the speed of light 

defined by the metric gab is unity. When u" is a unit vector, the square of the twist 
^a = ^abcdu'^^'^u'^ is a Combination of the Ci, C3, and C4 terms in the action ([3]), 

ujauj'' = -(V,Mb)(VV) + (V,M,)(VV) + (m^V.mJKVX). (6) 

In spherical symmetry the aether is hypersurface orthogonal, hence it has vanishing 
twist, so the C4 term can be absorbed by making the replacements 

ci -)■ ci + C4, C3 -> C3 - C4, C4 -)■ 0. (7) 

The field equations from varying (jS]) plus a matter action (coupled only to the 
metric) with respect to g*^^, u" and A are given by 

Gab = tI^^ + MGT'\b (8) 

gabU'^u' = 1, (10) 

where 

JV = K'^'mnVbU''. (11) 

The aether stress tensor is given by 

+ ci [(V^n,)(V"M,) - {VaUm){Vbun] 

+ C4 UaUb 

+ [Mn(VmJ"'") - C4M^] UaUb 

- ^Lugab, (12) 

where L^ = —K°'''mn^aU"^^bu"'- The Lagrange multiplier A has been eliminated from 
( !T2|) by solving for it via the contraction of the aether field equation (Q with m". 
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3.1. Metric redefinitions 

It is sometimes convenient to re-express the theory in terms of a new metric and aether 
field, related to the original fields by a field redefinition of the form 

9'ab = 9ab + (cr - l)^^^;, (13) 

u'" = -^u". (14) 

'a 



The constant a is restricted to be positive so the new metric remains Lorentzian. In 
effect, the field redefinition "stretches" the metric tensor in the aether direction by a 
factor 0". The action ([3]) for {g'ah, u'"') takes the same form as that for {gab, u"-) up to the 
values of the coefficients Cj. The Cj coefficients for the action expressed in terms of the 
new fields ([I3])-([T1]) are pTO] 

c'l = f ((1 + ^~')ci + (1 - <y-')c, - (1 - o-y) (15) 

4 = a(c2 + 1 - a-1) (16) 



„/ ^ 



C3 = 2 ((1 - ^"')ci + (1 + ^"')^3 - (1 - <y-')) (17) 

4 = C4 - |((1 - a-^)% + (1 - a-2)c3 - (1 - a-^f). (18) 

Certain combinations of the coefficients change by a simple scaling under the field 
redefinitions: 

Ci4 = Ci4 (19) 

C'l23 = ^Ci23 (20) 

C'lg - 1 = a(ci3 - 1) (21) 

c;-c'3-l=(7-l(Ci-C3-l). (22) 

The field redefinition relates solutions to the field equations coming from the two actions. 
Using a field redefinition the general form of the action can be simplified [T6] by 
eliminating one of the Cj or a combination of the Cj. If we choose a = (^2)^ = 1/(1 — C13), 
then c'^3 vanishes, i.e. C3 = —c[. In this case the two corresponding terms in the 
Lagrangian combine to make a Maxwell-like lagrangian, for which the Levi-Civita 
connection drops out. (This choice of a is positive and therefore preserves Lorentzian 
signature provided the original coefficients satisfy C13 < 1.) In the context of spherical 
symmetry, we may also exploit the vanishing to the twist of the aether to absorb the C4 
term in (|3]) by the replacements ([7]) as explained above. After these two changes, the 
Lagrangian takes a much simpler reduced form characterized by only two Cj coefficients, 

L.e = ~R- ^F-'Fab - c'^iVaU'^f - XiQabU'^u' - 1), (23) 

where Fab = 2V[aMb]. In the next section we will use this reduced form of the 
action to investigate the general behavior of stationary spherically symmetric solutions, 
addressing the existence of solutions around a metric horizon, asymptotic fiatness, and 
the regularity of the spin-0 horizon. 
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Another useful choice of field redefinition is to arrange for the new metric in fll3p to 
coincide with the effective metric for one of the wave modes ([2]) by choosing a = sj. In 
this way we transform to a "frame" where one of the spin-2, spin-1, or spin-0 horizons 
coincides with the metric horizon. Under the field redefinition all the squared speeds 
become s'^^ = sf/a, since the metric tensor in the aether direction is stretched by 
the factor a. In section [5] we use this method to make the spin-0 and metric horizons 
coincide, which simplifies the expansion of the field equations around the spin-0 horizon. 

4. Generic behavior of horizons and spatial infinity 

In this section we demonstrate that, for generic values of the q coefficients (at least 
for the reduced theory ( l23l) ). (i) metric horizons are generically regular in the three- 
parameter family of local solutions, (ii) asymptotic flatness imposes one condition on 
this family, and (iii) regularity of the spin-0 horizon imposes another condition, leaving 
a one-parameter family of black hole solutions with regular spin-0 horizons. 

4.I. Metric horizon expansion 

One way to determine the number of independent solutions with a regular metric horizon 
is to expand the field equations in a power series about a candidate horizon and solve 
algebraically order by order. We did this using the Maple computer application to carry 
out the algebra. The computation was prohibitively complicated using our methods 
with general values of the Cj, so we restricted attention to the reduced theory (l23l) . 
The determination of the solution space can also be done more "experimentally," by 
numerical integration of the field equations with varying initial data. With the latter 
method no special restriction on the c, is required. 

For the study of spherical black holes we adopt Eddington-Finkelstein (EF) type 
coordinates {v, r, 6, ip) with line element 

ds"^ = N{r)dv^ - 2B{r-)dvdr - r'^dVt^ (24) 

with advanced (null) time coordinate v and radial "area coordinate" r. The time- 
translation Killing vector is d^, and a metric horizon corresponds to N{r) = 0. These 
coordinates are regular at metric and other horizons so are useful for studying black 
holes and their interiors. Using these coordinates a stationary spherical aether field 
takes the form 

m" = a{r)d^, + b{r)dr (25) 

and the unit constraint (1101) becomes 

Na'^ - 2Bah = 1. (26) 

The field equations ([8]) and Qj become a set of coupled, second order ordinary differential 
equations (ODE's) involving the functions A^, B, a, b. We use the constraint (l26l) to solve 
for b in terms of the other three functions. Even for the reduced case C13 = C4 = the 
equations are sufficiently complicated that it does not seem useful to display them here. 
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Regularity at the metric horizon can be imposed by making a power series expansion 
about the radius r^ where A^(r^) = 0, 

N{r) = N'{rh){r - r^) + lN"{rh){r - rj,)^ + ■■■ (27) 

B{r) = B{rh) + B'{rh){r - m) + \B"{rh){r - r^f + ■■■ (28) 

a{r) = a{rh) + a\rh)ir-rh) + ^a"{rh)ir -Thf ^ (29) 

Inserting these expansions in the field equations, one can solve order by order for the 
power series coefficients. This allows the set of free parameters in the initial data at 
the horizon to be identified. At zeroth order in (r — Vh) the field equations imply that 
a'(r/j) is a function of N'{rh), B{rh), a{rh), and r^. The specific result is sufficiently 
complicated that it too does not seem useful to display here. Solving to higher orders we 
find that all remaining coefficients in the series expansion are determined by these four 
initial data parameters. Using the scaling freedom in the v coordinate {v — )• Xv where A 
is a constant) one of the initial values at the horizon can be fixed arbitrarily. Thus, there 
is a three-parameter family of local solutions with a regular metric horizon. As discussed 
in [13], we also found a three-parameter family of local solutions expanding about an 
arbitrary radius (i.e. not imposing a horizon), hence we conclude that regularity of the 
metric horizon generically imposes no restriction on the solutions. 

4-2. Asymptotic expansion 

In addition to regularity at all the horizons, we require the black hole solutions to be 
asymptotically fiat. To determine the form of such solutions we can change to the inverse 
radius variable x = 1/r and expand around x = 0, as was done in [1] (where isotropic 
coordinates were employed). In the reduced theory (!23|) this yields the solutions 

N{x) = 1 + Nix + —c[N^x^ + ■■■ (30) 

4o 

Bix) = 1 + ^c[N^x' - ^c[N',x^ + ■■■ (31) 

lb 1/ 

a{x) = 1 - ^Nix + a2X^ + {-—c[Nf + —Nf - N^a2)x'^ + ■■■ (32) 

9o lo 

where A^^i = N'{x = 0) and 02 = a"{x = 0), and the freedom to rescale v has been 

exploited to set N{x = 0) = 1. No more free parameters appear at higher orders, so the 

asymptotically fiat solutions are determined by the two free parameters A'^i and 02- 

An asymptotically fiat solution can be determined by using a simple shooting 

method, numerically integrating outward from an interior radius where there are three 

free initial data parameters. As in [13], we find that to match the asymptotic form 

( I30|) -(l32l) requires tuning just one of the three initial parameters, as expected since 

there are two free parameters in the asymptotic form. We conclude that in particular 

there is a two-parameter family of asymptotically flat "black hole" solutions with metric 

horizon fixed to lie at a given radius Vh- In practice, to integrate outward from a metric 

horizon we found it necessary to first use the perturbative solution about the horizon 
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to generate from the horizon data an initial data set some small radial distance away. 
This is because the ODE's have a singular point at the horizon. 

A more direct way to generate such asymptotically flat black hole solutions is to 
start the numerical integration near infinity and integrate inward using the inverse radius 
coordinate x. Since x = is a singular point of the ODE, it is necessary to start the 
integration at some small non-zero x value. The expansions f p0|) - fl32|) can be used to 
generate initial data as a function of Ni and 02- We find again this way that regularity 
at the metric horizon does not impose any conditions on A^i and 02- The functions A^, 
B, and a evolve smoothly through a point where A^ goes to zero. 

4-3. Asymptotically flat solutions and spin-0 horizon regularity 

So far we have shown that there is a two-parameter family of asymptotically flat black 
hole solutions with a regular metric horizon. Normally we expect just one black hole 
parameter, the total mass, unless there are conserved charges that can be additional 
parameters. In ae-theory there seems to be no such conserved charge, so the situation 
is puzzling. Another puzzling aspect of these black hole solutions not yet discussed here 
is that some have internal singularities at nonzero radius, rather than just at r = like 
most known black holes. Moreover, in some solutions we found that such singularities 
can occur externally, i.e. not inside a metric horizon. All these puzzles are resolved by 
the recognition that the singularities in question occur precisely at the location of the 
spin-0 horizon. Imposing regularity at the spin-0 horizon eliminates one free parameter, 
leaving us with a conventional one-parameter family of asymptotically flat black holes. 

In the rest of this subsection the full range of behavior of asymptotically flat 
solutions is discussed. In particular it is demonstrated that when a spin-0 horizon 
occurs it is singular for generic values of the two initial data parameters A^i and 02 at 
infinity. Evidence is given that by tuning one of these parameters to a special value a 
regular spin-0 horizon can be obtained. In the next section we show by a power series 
expansion around the spin-0 horizon that such regular solutions do indeed exist. 

For the reduced theory ( 123|) . the spin-2 and spin-1 speeds are both unity, but the 
spin-0 mode has squared speed Sq = {c2/c[){2 — c[)/{2 + 802) (relative to the aether), 
which is generically different from unity. At the spin-0 horizon the mode propagates at 
fixed radius, hence the surface of constant r at that location is null with respect to the 
effective metric g^i^ defined in (J2]). This implies the condition g^^ = in EF coordinates, 
which we find occurs where 

Na = or . (33) 

1 + So I- So 

The first root is less than 1 so according to fl2^ occurs when 6 < 0, i.e. when the aether 
tips inward. For the second root the aether tips outward. The combination / = Na"^ 
is independent of the arbitrary scale for the v coordinate and equal to one at infinity. 
Inserting the expansions (l30l) and (l32l) we find the asymptotic form 



fix) = N{x)a{xy = 1 + 2{a2 - In^)x^ - Ni{a2 - In^)x^ + 0{x^) (34) 
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Curiously, this expansion is independent of both c[ and Cj through order x^ (but not 
beyond) and depends hnearly on 02 — ^Nf through order x^. This pattern suggests an 
analytic solution may be possible, but we will not pursue this here. The static aether 
solution studied in [T3] corresponds to the case where f{x) = 1 for all x, which occurs 
when aa/A^f = 3/8. 

As in the previous subsection, we study solutions obtained by integration inwards 
starting from an asymptotically flat spatial infinity. Since the theory has no length 
scale, the solution with data (A''i,a2) is trivially related to that with data (AA''i,A^a2) 
(as are Schwarzschild solutions with different mass trivially related). If we think of the 
line element ds"^ as giving a numerical value specified with respect to a given length 
unit, then to go from one solution to another we need only change the unit of length. 
Thus without loss of generality we can fix units with A'^i = ±1. The solutions then 
depend on the choice of theory through c'^, c'2, on the parameter 02, and on the sign of 
A^i. A systematic study of these solutions is beyond the scope of this paper; here we 
just indicate the various behaviors we have encountered, and then focus on the regular 
positive mass black holes. 

Let us first consider positive mass solutions, i.e. A^^i = — 1. As the radial coordinate 
decreases, A^ decreases from 1, while the combination Na^ increases or decreases 
according as 02 is greater or less than 3/8. There are solutions where /(r) does not 
reach or (1 =1= so)/(l =t -So) and therefore neither a metric nor spin-0 horizon is attained. 
In some of these N{r) re-curves out to positive infinity, a{r) approaches zero and the 
solution reaches a curvature singularity near r = 0. There are also solutions similar to 
the static aether of [13], where B{r) goes to infinity as N{r) approaches a finite value, 
indicating a minimal area two-sphere. In some cases the larger root for a spin zero 
horizon in fl33|) is reached by f{r). In other solutions /(r) may reach a metric horizon, 
but does not attain the value corresponding to a spin-0 horizon. This can only happen 
when So > 1. N{r) again re-curves out to infinity and a{r) approaches zero near r = 
and there are outer and inner metric horizons. In still other solutions, the functions and 
their derivatives are regular up to the point where /(r) reaches the spin-0 horizon, but 
generically the spacetime is singular at that point. If Sq > 1, /(r/j) is negative and the 
singularity is located inside a metric horizon. If sq < 1 the singularity occurs without 
any metric horizon. 

For a specific example of this last type we choose parameters c[ = 0.051 and 
C2 = 0.116, for which the spin-0 speed is 1.37. (These arise from starting with coefficients 
that satisfy all the observational constraints described in [6] and performing the field 
redefinition to the reduced action ( !23l) .) Figure [1] shows the behavior of A^, B, and a 
for A^^i = — 1 and 02 = —0.1. In this case there are "outer" and "inner" metric horizons 
where A^ = 0, but /(r) does not reach as low as —0.158, which is required by (133!) in 
this case for a spin-0 horizon. For values of 02 < —0.1 the minimum value of A^ shifts 
upward and eventually A^ never reaches zero, i.e. the metric horizon disappears. On the 
other hand, for 02 > —0.1 the minimum value of / decreases until the spin-0 horizon is 
reached. At this point a' goes to negative infinity, A^' blows up to positive infinity, and 
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Figure 1. Solution for reduced theory ([23|) with c'l — 0.051 and c'2 = 0.116, determined 
by data A^i = — 1 and a2 = —0.1 at spatial infinity. There is both an outer and inner 
metric horizon where N vanishes, but / does not decrease enough to reach a spin-0 
horizon. The functions N and / go to zero slightly inside r = 1 which would be 
the horizon radius of the corresponding Schwarzschild solution. As 02 increases the 
minimum of f{r) decreases until the solution acquires a spin-0 horizon, where is it 
generically singular. 



there is a curvature singularity. In contrast, note that in Figure [T] a' (r) has a maximum 
value while N'{r) goes to negative infinity. This suggests that at some special value of 
02 there is a transition were the concavity of a'{r) and N'{r) changes and the derivatives 
are finite at a spin-0 horizon. Regularity at the spin-0 horizon seems thus to impose one 
condition on the asymptotic values A^^i and 02- 

In the negative mass case one might expect only solutions analogous to negative 
mass Schwarzschild, with N increasing from 1 at spatial infinity and no spin-0 or metric 
horizon. While the solution does take this form for all 02 in the theory with c'l = 0.051, 
C2 = 0.116, for other values of c'12 there are ranges of 02 where N increases from 1 at 
infinity, but then decreases to a metric horizon at finite r, and all the functions and 
their derivatives are regular until /(r) reaches the value less than zero required for a 
spin-0 horizon. This peculiar behavior of a negative mass solution with metric and 
spin-0 horizons remains to be studied more closely. In particular, it is not clear whether 
a negative mass solution with a regular spin-0 horizon could exist. 



5. Black holes with regular spin-0 horizons 

In this section we discuss the behavior of black hole solutions possessing regular spin-0 
horizons. Rather than imposing regularity at the spin-0 horizon by the shooting method 



C2 = r.:\ .T:\. ' • (35) 
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integrating in from infinity, we instead expand the field equations in a power series about 
a non-singular spin-0 horizon. 

5.1. Horizon expansion 

Due to the complexity of the field equations and their singular nature at the horizon, 
we were unable to implement the power series solution about a spin-0 horizon in the 
generic reduced theory ( !23ll (even with computer aided algebra). It might be possible 
to obtain the perturbative solution by a more well-adapted method, but instead we 
simplified the computation by making a field redefinition to a new metric for which the 
spin-0 and metric horizons coincide. Starting from an arbitrary set of coefficients Cj, 
this is implemented by the choice a = Sq in (IHj), after which we have sq = 1 without 
loss of generality in the theory. As before we can also then exploit spherical symmetry 
to absorb C4 by making the replacements ([7]), which do not disrupt the coincidence of 
the spin-0 and metric horizons since this is just a re-expression of the same Lagrangian 
without changing the field variables |s|. This reduces the distinct parameter space to 
just (ci,C3) (omitting the prime in the notation for Ci^s). After this field redefinition the 
coefficient C2 is given by 

— ^C3 — C^ — ZO^C-^ — '-'I'-'S 

2-4ci + 3cf + ScgCi 

A further simplification of the equations is achieved by trading the metric function 
A^ for the combination of metric and aether functions / = Na"^. We have no insight 
into why this simplifies the expansion of the field equations about the common metric 
and spin-0 horizon at A^ = = /, although as stated above the combination Na'^ is 
invariant under a rescaling of the v coordinate. The field equations in this set of field 
variables involve a, a', a", /, /', /", B, and B'. At the horizon /(r) vanishes linearly, 
/(^) = f'{f'o){f' ~ '''o) + ■ ■ ■ By a constant rescaling of v we can furthermore set B{rh) 
equal to 1. Using this along with ( l28l) and ( l29ll the field equations can be expanded and 
solved order by order for the coefficients of the power series. 

Solving the field equations for this theory as algebraic equations for the expansion 
coefficients we find that at zeroth order in (r — Vh) the quantities a{rh), a'{rh), a"{rh), 
and f"{rh) are determined by free parameters r^, f'ijh), B'{rh). We succeeded in solving 
the equations to the next order in (r — r^) only in the special cases C3 = 0, C3 = ci, and 
C3 = — Ci. In these cases we find that B'{rh) is determined by rh and f'{rh). Hence, 
consistent with the expectation of the previous section, there is a two-parameter family 
of local solutions around the regular spin-0 horizon. These solutions are generically not 
asymptotically flat. 

S The spin-0 speed is invariant under ([7|), as guaranteed by this argument. The spin-1 speed is not 
invariant, but this does not contradict the argument since there is no spherically symmetric spin-1 
mode. Note however that in diagnosing whether spin-1 perturbations are trapped in a given black hole 
it is important to use the value of the spin-1 speed written in ([l} before the C4 coefhcient has been 
absorbed. 
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Figure 2. Plots of /, N, B, a, and 5* (the Schwarzschild version of N) vs. z = r ~ rh 
for ci = 0.3, in units with rp = 2. The horizon radius is r^ ~ 2.07 for the ae-theory 
black hole, so S = Q a,i z ~ —0.07. The solutions agree closely outside the horizon. 
Deviations are noticeable near the horizon and become significant in the interior, where 
N blows up more rapidly. Near the singularity / begins to oscillate rapidly. 



5.2. Asymptotically flat black holes 

To produce asymptotically fiat solutions we numerically integrate outward, starting with 
the horizon data and matching onto (ISUl) . fl?T]) . and fl52]) by tuning /'(r/^) until /(r) is 
constant and equal to 1 at very large r values. The asymptotic flatness boundary 
condition at inflnity thus reduces the number of free parameters to one, namely the 
horizon radius itself. Solutions with different horizon radii are trivially related. Since 
rh is a singular point of the ODE's, it is necessary to start the integration with initial 
data at some small positive value oi r ~ r^. We used the series solution determined by 
a given rh and f'{rh) to generate this initial data. To examine the solution inside the 
horizon, we numerically integrated inward, starting at a small negative value oi r — rh 
with data generated by the same series solution. 

Here we will discuss the properties of the solutions for the C3 = theory only, whose 
behavior is typical of the three special cases C3 = 0, ±ci. Figure |2] displays the solution 
for ci = 0.3, C2 = —.025, 03 = 04 = 0, together with S{r) = 1 — 2/r, the Schwarzschild 
version of N{r) with the same mass. For this plot we use the scaling freedom of v to 
convert the numerical solution to a "gauge" where the metric functions and v component 
of the aether are all equal to 1 at inflnity. The two metric functions B{r) and N{r) 
in GR and ae-theory are in very close agreement outside the horizon, while inside they 
differ noticeably. 
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5.2.1. Black hole mass The ADM mass M^^m of an asymptotically flat spacetime 
whose asymptotic metric takes the Schwarzschild form at 0(l/r) is directly determined 
by the coefficient vq = 2GMadm of the 0(l/r) part of gu- In ae-theory the relation 
between Madm and the total energy E of the spacetime is GMadm = GnE [171 [9], 
where Gn = G/{1 — C14/2) is the Newton constant appearing in the force law between 
two weakly gravitating masses [12]. We shall refer to the quantity ro/2 with dimensions 
of length as the "mass" in what follows, and denote it by M. For a Schwarzschild black 
hole in GR, tq is equal to the horizon radius r^. In ae-theory the ratio ro/r/^ is a constant 
(since there is only one length scale) determined by the coupling coefficients Cj. 
The EF line element (1211) transforms to Schwarzschild form 



ds^ = Ndt^- {B^/N) dr^ - r^dn^ (36) 

with time coordinate t defined by dt = dv — [B/N)dr. The asymptotic form of B 
( 13T|) shows that B = 1 + 0(l/r^), so up through 0(l/r) the line element (l36l) has the 
standard asymptotic form if A^ and B are normalized to 1 at infinity. In generating an 
asymptotically fiat numerical solution we fixed the scale freedom of the v coordinate 
by imposing B{rh) = 1 at the horizon however, so the asymptotic form of A^ is 
A^oo + Ni/r + 0(l/r^). The mass is given by M = ro/2 = Ni/2Noo, which can be 
extracted from the numerical solution at large r. 

5.2.2. Horizons The solution displayed in Figure 121 has metric and spin-0 horizons at 
2; = 0, but how about spin-1 and spin-2 horizons? Is the fastest speed actually trapped? 
The condition for a horizon corresponding to a speed sq is given in ( l33l) . As the speed 
approaches infinity the horizon value of / approaches —1 from above. In Figure 121 (and 
for all values of ci that we studied up to 0.7), the minimum value of /(r) is less than 
—1, which is sufficient to trap any wave mode. The fact that /(r) curves back to being 
greater than —1 indicates that an inner horizon might exist for some wave modes in 
certain parameter ranges of q. 

In the theory under discussion we have C3 = C4 = and C2 = —c\/(2 — 4ci + 3c^), so 
the squared mode speeds in ^ are given by 1/(1 — ci) for spin-2 and (1 — ci/2)/(l — ci) 
for spin-1. With < Ci < 1 both of these are greater than 1, and the spin-2 speed 
is the highest. In the particular case shown in the figure, the spin-1 speed is 1.10 and 
the spin-2 speed is 1.20, which correspond to horizons at / = —0.049 and / = —0.089 
respectively, which do not seem to be reached a second time. 

5.2.3. Oscillations A notable aspect of the black hole interior displayed in Figure is 
the oscillation in /. The function h = Ba also oscillates in a similar manner, but is 180 
degrees out of phase. In addition, there are related oscillations in the curvature scalar 
and aether congruence behavior discussed below. These oscillations are reminiscent 
of the interior behavior found in Einstein- Yang- Mills black holes [18], where the metric 
functions and derivative of the Yang- Mills potential oscillate an infinite number of times 
before the singularity. 
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Figure 3. Oscillations of / = Na^ near the singularity inside the black hole solution 
shown in figure [21 plotted vs. z = r — r^ in units with A/ = 1. 



While N decreases monotonically and B increases monotonically, a goes to zero, 
so the oscillations of Na^ and Ba arise because of variations in the magnitude of their 
derivatives. Since the oscillations inside z = —2 are not clearly visible in Figure [2] a 
zoomed in graph of f{z) is provided in Figure [31 From this graph it is clear that / 
smoothly turns over at least once more before the singularity. Although the number 
of oscillations before the singularity appears finite, it is possible that the numerical 
integration employed is not capable of resolving additional or even infinitely more 
oscillations. More information may be obtained in the future by improved numeric 
methods or analytic methods around r = 0. 



5.2.4- Curvature singularity There appears to be a spacelike curvature singularity at 
or near r = 0, as in the Schwarzschild solution of GR. In Figure [2, the approach of 
A^ to negative infinity near r = suggests a singularity. In Figure [H the logarithm of 
the Kretschmann scalar K = Rabcd.R"'^'^'^ is plotted vs. Inr for the ae-theory solution 
together with its value in the corresponding Schwarzschild solution with the same mass. 
In the latter case K = 48/r^ in units with M = 1, so logK = —6 Inr + In 48. The rate 
dlnK/d\nr for the ae-theory solution seems to alternate between roughly —6 and —4.5. 
The location of the transitions may be correlated with the oscillations discussed above. 
The location of the singularity seems to be at r = for all the values of ci, although 
our numerical solutions do not permit a determination of the exact location. 
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Figure 4. Plot of \n RabcdR vs. Inr for ci — 0.3 ae-theory black hole (wiggly 
curve) and Schwarzschild black hole (stright line) of the same mass, in units with 
M = 1. r is plotted on a logarithmic scale. The slope for the GR case is —6, while for 
the ae-theory case it alternates between roughly —6 and —4.5. 



5. 2. 5. Aether congruence The aether field defines a congruence of radial timelike curves 
at rest at infinity and fiowing into the black hole. It is interesting to compare this with 
the static frame and with the congruence of freely falling radial geodesies with 4- velocity 
v" that are also at rest at infinity. Being unit vector fields, at each point u"" and f " can 
be fully characterized by their Killing energy, i.e. their inner product with the Killing 
vector. The free-fall congruence has a conserved energy that is equal to one if the Killing 
vector is normalized to one at infinity. The aether does not fall as quickly outside the 
black hole. In fact it remains rather aligned with the Killing vector up until quite close 
to the horizon. 

To characterize and contrast the free-fall and aether congruences we plot in Figure 
|5]the derivative dr/dr of radius with respect to proper time along each congruence. Let 
us call this the quantity the "proper velocity" IjjJ The aether and free-fall are both at 
rest at infinity, but only as the horizon is approached is the aether finally pulled away 
from the Killing direction. As close as r = Sr^ {z ^ 4), the proper velocity of the aether 

II The magnitude of this quantity is affected both by the radial motion and the behavior of the proper 
time. For instance as the particle becomes lightlike the proper time goes to zero and this derivative 
diverges. However we could think of no better measure of the radial velocity. One might use the 
3-velocity relative to a static observer outside the black hole, but since the static observer becomes 
lightlike at the horizon, this 3-velocity will be equal to one at the horizon for any finite timelike 4- 
velocity, so it does not distinguish different motions there. 
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Figure 5. Radial proper velocity dr/dr of free-fall (lower, solid curve) and aether 
(upper, dashed curve) vs. z — r ~ Vh, in units with M = 1, for ci = 0.3. In contrast 
to the free-falling geodesies, the aether does not begin to fall significantly inward until 
close to the horizon. 



is still about fifteen times smaller than that of free-fall. Inside the horizon the aether 
proper velocity is equal to the free-fall one around z = —1.3, but the 4-velocities do 
not agree there. The aether is still going inward faster, but its proper time is "running 
slower" so it can have the same proper velocity. 

To compare the aether and free-fall motions inside the horizon we plot in Figure [6] 
the inward 3-velocity of the aether with respect to the free-falling frame. The relative 
velocity is initially negative, meaning that the aether is not falling in as fast as the the 
free-fall frame. It is clear from this plot at around z = —1.3 the aether still lags well 
behind free-fall. However, around z = —1.9 the relative velocity is zero, and after that 
it oscillates a couple of times (at least) before reaching the singularity. 



5.2.6. Surface gravity and the first law of black hole mechanics The laws of black hole 
mechanics have been shown to apply to a wide class of generally covariant metric theories 
of gravity coupled to matter [19] . There appears to be no straightforward extension of 
the first law and the concept of black hole entropy to ae-theory however [17j, a difficulty 
that is tied to the fact that there is no smooth extension of the aether to the bifurcation 
surface of the Killing horizon. Moreover, it is not clear to which horizon the law should 
apply, in a theory with multiple characteristic surfaces. For example, in the solutions 
considered in this section, the spin-2 horizon is inside the spin-1 horizon which is inside 
the joint spin-0 and metric horizon. One might imagine that the relevant horizon is 



Black Holes in Einstein- Aether Theory 

0.4- 



19 




-0.5 



Figure 6. Inward 3-velocity of the aether relative to free-fall inside the horizon for 
ci — 0.3. The velocity is initially negative and the aether lags behind the free-fall. 
Near the singularity the velocity oscillates between between faster and slower than 
free-fall. 



always the Killing horizon, but recall that by a field redefinition we can make any one 
of these horizons be the Killing horizon. We shall not try to shed any light on these 
puzzling issues here. Rather, we just briefiy examine the variational relation between 
mass, surface gravity and area of the spin-0 horizon, for possible future use. 

The first law of black hole mechanics for spherically symmetric neutral black holes 
in GR takes the form 

kSA 



6M 



a 



SttG" 



(37) 



where A = 47rr^ is the horizon area, k is the surface gravity, and a = 1. By dimensional 
analysis such a variational relation must also hold in ae-theory, with some value for 
the dimensionless constant a that depends on the dimensionless coupling coefficients 
Cj. Presumably for M we should put the total energy E of the spacetime, and for "G" 
we should put the Newton constant Gn governing the attractive force between distant 
bodies. Alternatively one might use the ADM mass Madm and the constant G appearing 
in the ae-theory action (J3]). As discussed in section 15.2.11 GMadm = GnE = ro/2, so 
these two choices actually yield identical "first laws" . If we express the mass and area 
in terms of tq and r^ respectively, ( 137|) thus becomes 

6ro = 2aKrhSrh. (38) 

As pointed out in section 15.2. H ro and Vh are proportional, so we infer that a is 
determined by the dimensionless combination 

(39) 



a 



2nrr 
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Table 1. Properties of black hole solutions for several ci values, in units with r^ — 1. 



Cl 


f'irh) 


7 = U"-Va 


''o 


K 


a 


0.1 


2.096 


1.619 


0.990 


0.507 


0.976 


0.2 


2.072 


1.608 


0.979 


0.517 


0.947 


0.3 


2.039 


1.592 


0.966 


0.528 


0.914 


0.4 


1.997 


1.568 


0.951 


0.543 


0.876 


0.5 


1.941 


1.535 


0.933 


0.562 


0.830 


0.6 


1.867 


1.490 


0.911 


0.588 


0.787 


0.7 


1.767 


1.429 


0.881 


0.625 


0.704 



which depends on the coefficients Cj defining the theory. 

5.2.7. Black hole properties for different values of Ci Various properties of the black 
hole solutions for different values of Ci are displayed in Table [H The other coupling 
coefficients have the values C3 = C4 = and C2 is given by (!35|) . For each ci there is 
a one-parameter family of black hole solutions with regular spin-0 horizon, labeled by 
horizon radius. For the values in the table we compare black holes with the same horizon 
radius, and adopt units with r/j = 1. The Killing vector which enters the definition of k 
and a is normalized to unity at spatial infinity. 

The values of f'{rh) that yield asymptotically fiat solutions for different choices of 
Cl are displayed in the 2nd column. These values decrease as Ci grows. For ci = 0.8 
and larger we could not find a fir^) that yielded an asymptotically fiat solution. The 
third column shows the gamma factor between the aether and free-fall velocity at the 
horizon. The fourth column shows Tq = 2GMj\dm- This is equal to r^ for Ci = (a 
Schwarzschild black hole), and decreases by 12% as q increases up to 0.7. Conversely, 
for a given mass the black hole horizon is larger for larger Ci. The fifth column shows 
the surface gravity, which for ci = is l/2r/i and increases by 25% as ci increases up 
to 0.7. The last column is the dimensionless ratio ( l39l) appearing in the first law ( 1371) . 
which is unity for ci = and decreases by 30% as Ci increases to 0.7. 

6. Discussion 

In this paper we considered the meaning of a black hole in Einstein-Aether theory, 
arguing that the fastest wave mode must be trapped if the configuration is to qualify as 
a causal black hole. Regularity at the spin-0 horizon was identified as a key property of 
black holes in ae-theory. It was found that, for generic values of the coupling constants 
Cj, regularity at a metric horizon imposes no restrictions on spherically symmetric, static 
local solutions but regularity at a spin-0 horizon imposes one condition. At least for a 
class of coupling constants, there is a one parameter family of asymptotically fiat black 
hole solutions with all horizons (metric and spin-0, 1,2) regular. 

From an astrophysical point of view an essential question is what happens when 
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matter collapses. It is a plausible conjecture that nonsingular spherically symmetric 
initial data will evolve to one of the regular black holes whose existence has been 
demonstrated here, but this has certainly not been shown. It would be very interesting 
to answer this question by numerical evolution of the time-dependent field equations. 
To do so, one could add scalar matter to form a collapsing pulse, but this is likely not 
necessary since the aether itself has a spherically symmetric radial tilt mode that can 
serve the purpose. This would correspond to formation of a black hole by an imploding 
spherical aether wave. 

We examined several properties of regular black hole solutions for a special class of 
coupling coefficients defining the ae-theory. A complete classification of the solutions for 
different coefficients remains an open research problem, as does the study of non-black 
hole solutions and negative mass solutions, which strangely enough could include black 
holes. Oscillating behavior approaching the internal singularity has been identified, but 
not studied in detail. Finally, only the static, spherically symmetric case was examined; 
the question of rotating black hole solutions remains untouched. 
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